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VM 型 脉冲 管制 冷 机 数值 模拟 与 实验 研究 
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摘 要 本 文 介绍 了 由 VM 型 热 压缩 机 驱动 的 脉冲 管制 冷 机 基本 结构 、 数 值 模拟 与 实验 研究 . 热 压缩 机 借鉴 VM 制冷 
机 产生 压力 波动 的 方式 , 即 依靠 室温 与 冷 源 之 间 的 温度 差 产生 压力 波动 ,驱动 低温 级 脉冲 管制 冷 机 . 利用 液 氮 温 区 制冷 机 
为 热 压缩 机 提供 所 需要 的 冷 量 ， 避 免 使 用 低温 液体 带 来 的 不 便 . 低温 级 脉冲 管制 冷 机 的 回 热 属 与 脉 管 采用 U 型 布置 , 调 
相机 构 采 用 小 孔 气 库 加 双向 进 气 方式 . 这 种 结构 不 使 用 有 降 压缩 机 ， 可 以 产生 低频 压力 波动 ,保证 回 热 属 在 低温 区 高 效 工 
作 . 在 原 有 结构 的 基础 上 , 改造 热 压缩 机 回 热 属 斥 十 与 填料 , 使 用 Sage 软件 进行 了 整 机 模拟 , 对 调 相 机 构 进行 优化 , 考察 
其 在 不 同 的 运行 频率 、 平均 压力 、 排 出 器 位 移 下 的 制冷 性 能 . 对 改造 后 的 脉冲 管制 冷 机 进行 了 初步 的 实验 研究 ,获得 14.4 
K 无 负荷 制冷 温度 . 


关键 词 VM; 热 压缩 机 ; 小 孔 ; 双向 进 气 ; 数值 模拟 .给 
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The Numerical Simulation and Experim Btuqy on a VM Type Pulse Tube 
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Abstract This paper mai 
study on a VM type th 
that in VM cryocoolerNuti Zes temperature difference between room temperature and cold heat Sink 


Vantroduces the structure, numerical Simulation and experimental 
-driven pulse tube cryocooler. The thermal compressor, similar to 
to generate pressure wave. Two cryocoolers working at liquid nitr temperature are thermally 
coupled with the cold end of this thermal compressor to pro *the required precooling power, 
preventing the inconvenience from cryogenic fluid. Low temperat'ite stage regenerator and pulse tube 
are arranged with U type structure, and the orifice andxt double-inlet are used as phase shifters. 


This cryocooler consists of a thermal compressor Wy of valved compressor and operates with a 


low frequency, which is beneficial to improving re rator performance at low temperature. Based 
on the previous study, the regenerator for thexifar compressor is modified, and numerical simulations 
on the pulse tube cryocooler are perfor The dependence of system performance on frequency, 
average pressure and displacer displacement amplitude are investigated in the simulation. Finally, pre- 


liminary experimental study were performed and a lowest no-load temperature of 14.4 K were achieved. 


Key words VM; thermal compressor; orifice; double-inlet; numerical simulation 
0 引 言 


随 着 科技 事业 的 不 断 发 展 ， 低 温 技术 占据 的 地 
位 愈加 关键 .小 型 低温 制冷 机 是 低温 技术 的 重要 分 
支 之 一 , 在 低温 电子 、 红 外 器 冷却 、 超 导 滤 波 、 核 磁 
成 像 、 强 电 超 导 等 场合 具有 重要 应 用 . 目前 使 用 比较 
广泛 的 小 型 低温 制冷 机 主要 有 : Joule-Thomson(J-T) 
制冷 机 、Gifford-McMahon(GM) 或 GM 型 脉冲 管制 
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冷 机 、Stirling 或 Stirling 型 脉冲 管制 冷 机 . JT 节 流 
制冷 机 冷 头 没有 运动 部 件 , 振动 小 , 可 靠 性 高 ; 但 由 
于 节 流 是 典型 的 不 可 道 过程， 内 课 效 率 低 是 ， 且 对 
工 质 洁净 度 要 求 非常 高 , G-M 或 GM 型 脉冲 管制 冷 
机 技术 发 展 成 熟 , 但 冷 端 振动 大 、 体积 大 、 需要 定期 
维护 ,采用 GM 压缩 机 与 配 气 阀 门 提供 压力 波 , 具 
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有 内 店 的 不 可 逆 损 失 . Stirling 型 制冷 机 基于 Stirling 
循环 , 具有 潜在 的 高 效 性 , 但 其 通常 工作 在 高 频 工 况 
下 , 回 热 器 在 低温 区 尤其 到 液 所 温 区 损失 很 大 , 效率 
较 低 . 为 了 获得 极 低温 区 , 需要 采用 三 级 及 三 级 以 上 
结构 ,系统 复杂 ,热效率 低 4. 

针对 不 同 应 用 领域 对 液 氨 低 温 制 冷 机 提出 的 更 
高 要 求 , 诸如 系统 紧凑 、 效 率 高 、 振 动 噪音 小 、 寿 合 
长 、 免 维护 周期 长 等 , 借鉴 VM 制冷 机 产生 压力 波动 
的 方式 回 ， 基 于 脉冲 管制 冷 机 冷 端 无 运动 部 件 、 结 
构 简单 、 可 靠 性 高 、 寿 命 长 等 优点 , 文献 [6] 发 展 了 
干 式 (cryogen-free)VM 型 热 压缩 机 驱动 的 脉冲 管制 
冷 机 .其 中 ， 热 压缩 机 的 排出 器 由 动 磁 式 直线 电机 
驱动 , 依靠 室温 与 冷 源 之 间 的 温度 差 产生 压力 波动 ; 
台 液 握 温 区 斯 林 型 脉冲 管制 冷 机 提供 热 压缩 机 所 
需要 的 冷 量 ， 有 利于 优化 控制 热 压缩 机 低温 端 温度 
和 实现 结构 上 的 紧凑 性 ;脉冲 管 制冷 机 由 热 压 缩 机 
产生 的 压力 波动 驱动 。 该 制冷 机 运行 频率 低 、 驱 动 
切 小 、 结 构 简单 ， 对 于 发 展 液 氨 温 区 高 效 紧 型 制冷 
机 具有 重要 的 研究 意义 

本 文 在 原 有 结构 的 基础 上 ,修改 热 压缩 
器 以 提高 其 输出 压 比 ,对 修改 后 的 制冷 
模拟 ， 考 察 其 在 不 同 频率 、 平 均 压 人 
移 下 的 性 能 ， semen. 


1 系统 结 :和 

在 文献 [7] 工作 的 基础 上 ， 对 热 压缩 机 回 热 器 
进行 改造 以 提高 其 输出 压 比 . VM 型 热 压缩 机 驱动 
脉冲 管制 冷 机 结构 如 图 1 所 示 ， 基 于 Sage 软件 的 
整 机 模拟 程序 如 图 2 所 示 , 图 2 左 半 部 分 为 热 压缩 
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机 模型 ， 右 半 部 分 为 制冷 机 模型 。 改 造 方案 基于 的 
思路 是 ; 先 单独 计算 低温 脉冲 管制 冷 机 ， 提 取 其 入 
口 阻抗 幅 值 与 相位 角 ; 将 其 阻抗 带 入 热 压 缩 机 程序 
匹配 计算 ; 以 效率 最 高 为 目标 函数 (效率 = 输出 声 
功 / 预 冷 量 )， 以 回 热 器 热 端 气体 位 移 幅 值 / 回 热 器 长 
度 和 20% 为 限制 条 件 , 优化 回 热 器 的 直径 与 长 度 ; 在 
最 优 直 径 与 长 度 下 ， 优 化 回 热 融 不 锈 钢丝 网 填料 的 
目 数 与 丝 径 . 经 过 计算 ， 最 终 选 定 其 参数 如 表 1 所 
示 ， 回 热 需 填 充 80 目 不 锈 钢丝 网 ， 丝 径 为 120 bm。 
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Fig. 2 Numerical model of VM type pulse tube cryocooler 
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表 1 热 压缩 机 运行 参数 及 回 热 器 结构 尺寸 机 回 热 右 热 端 、 冷 端的 声 功 随 平均 压力 的 变化 趋势 。 
Table 1 The operating parameters of thermal 其 中 ， 热 端 声 功 为 图 1 所 示 的 Eo,RegHot， 冷 端 声 功 
compressor and structure dimensions of 为 图 1 所 示 的 Fz,RegCold. 随 着 平均 压力 的 增 大 , 热 
Wee 压缩 机 压力 波动 幅 值 增 大 ， 输 出 到 回 热 器 热 、 冷 端 
Pn /MPa f/Hz D/mm L/mm 声 功 也 增 大 ， 
2.5 3 30 150 
04 
整 机 数值 模拟 | 
en 
在 确定 的 热 压 缩 机 回 热 器 结构 尺寸 上 ， 对 VM ee 
热 压缩 机 型 脉冲 管制 冷 机 进行 整 机 模拟 计算 ， 考察 a | 昌 
不 同 运行 频率 、 平 均 压力 与 热 压缩 机 排出 器 位 移 下 和 |“ 咒 
制冷 机 的 性 能 . 制冷 机 主要 部 件 的 尺寸 如 表 2 所 示 ， ae 人 
: : : 10.1 地 
表 2 主要 部 件 参数 a 
Table 2 Parameters of main components | 


14 16 18 20 22 24 2.6 


一 主要 部 件 参数 
> 制冷 机 D=15 mm, L=210 mm, HoCu2， 平均 压力 /MPa 
pe es RS 热 压 缩 机 压力 波动 幅 值 与 制冷 机 回 热 器 热 端 、 冷 端 声 功 
nn ig. 4 Dependency of pressure wave amplitude and acoustic 
ower at warm and cold end for low stage regenerator on the 
2.1 不 同 平均 压力 下 的 计算 average pressure 
该 计算 中 ,保持 频率 3 Hz， nd 
热 压 缩 机 冷 热 问 温度 分 别 为 300 K 和 器 下 改变 二 图 5 显示 的 是 热 压缩 机 冷 端 所 需 预 冷 量 随 平 


均 压力 ， 在 不 同 平均 压力 下 优化 湖 演 糙 必 向 进 气 网 。 均 压力 的 变化 站 势 ， 随 着 平均 压力 的 升 高 ， 所 需 机 
门 开 度 , 获得 最 低 制冷 温度 . 形 斤 角力 在 1.2 MPa 到 冷 量 不 断 增 大 ， 与 输出 声 功 随 平均 压力 的 变化 趋势 
2.5 MPa 之 间 变 化 . \ 

图 3 显示 的 是 制冷 机 最 低 制 冷 温度 随 平 均 压 力 
的 变化 趋势 。 随 着 平均 压力 的 增 大 ,无 负荷 制冷 温 
度 降低 . 考虑 到 He-4 物性 随 压力 急剧 变化 ， 以 及 系 
统 耐 压强 度 的 限制 ,其 平均 压力 不 能 过 高 


== 
CE 
[¢ 


a 


温度 /KK 


1.8 2.0 22 2.4 
平均 压力 /MPa 

图 5 热 压 缩 机 所 需 预 冷 量 随 平均 压力 的 变化 

| | | | | | | ( 冷 端 温度 固定 为 90 K) 

10 12 14 16 18 20 22 24 2.6 Fig. 5 Dependency of required precooling power on the 


平均 压力 /MPa average pressure 
图 3 最 低 制 冷 温度 随 平均 压力 的 变化 a i 
Fig. 3 Dependency of lowest no-load temperature on the 图 6 显示 的 是 制冷 机 回 热 器 与 脉 管 冷 、 热 端 阻 
average Pressure 抗 相位 角 随 平均 压力 的 变化 趋势 . 


位 关系 是 在 其 中 部 某 处 压力 波 与 质量 流 同 相 四. 从 
图 4 显示 的 是 热 压 缩 机 压力 波动 幅 值 以 及 制冷 a 
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系 , 主要 原因 在 于 肪 管内 部 相位 变化 跨度 大 ， 由 其 
空 体积 导致 ， 这 对 日 后 制冷 机 的 改进 与 优化 具有 一 
定 指导 意义 。 


加 可 热 ( 噶 热 问 端 


相位 角 /() 


Te 0 2 2 
平均 压力 /MPa 
图 6 制冷 机 回 热 器 、 脉 管 冷 热 端 阻抗 相位 随 平 均 压 力 的 变化 


Fig. 6 Dependency of phase difference between pressure wave 


and mass flow rate at warm and cold end for low stage 


regenerator and pulse tube on the average pressure 2 


2.2 不 同 频 率 下 的 计算 

制冷 机 平均 压力 2.5 MPa， 排 出 i 
gp 
7 获得 最 低 制 
冷 温度 ， 

0 
趋势 . 随 着 频率 的 升 高 , 无 负荷 温度 先 减 小 后 增 大 ， 
在 4 Hz 下 达到 最 低 . 从 图 8 发 现 , 热 压缩 机 压力 波 
动 幅 值 与 制冷 机 回 热 器 热 端 声 功 均 随 着 频率 的 升 高 
而 增 大 . 相 比 其 他 情况 ,1 Hz 下 由 于 输出 声 功 急剧 
减 小 , 无 负荷 温度 迅速 升 高. 


频率 /Hz 
图 7 最 低 制 冷 温 度 随 频率 的 变化 


Fig. 7 Dependency of lowest no-load temperature on the 


frequency 
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图 8 热 压缩 机 压力 波动 幅 值 与 制冷 机 回 热 器 热 端 、 冷 端 声 功 随 
频率 的 变化 


Fig. 8 Dependency of pressure wave amplitude and acoustic 
power at warm and cold end for low stage regenerator on the 
frequency 


AAA 

频率 从 1 Hz 变化 到 6 Hz, 输出 声 功 增 大 

后 时 ， 所 需 预 冷 量 从 8.12 W 变化 到 58.17 W. 由 
于 预 冷 量 是 由 相对 固定 输入 功率 的 脉 管制 冷 机 提供 ， 
所 需 预 冷 量 增 大 使 得 热 压缩 机 冷 端 温 度 升 高 ， 减 少 
了 其 工作 温差 , 对 性 能 有 一 定 负面 影响 . 
2.3 不 同 排出 器 位 移 下 的 计算 

制冷 机 平均 压力 2.5 MPa, 频率 3 Hz， 热 压缩 
机 冷 热 端 温度 分 别 为 300 K 和 90 K, 排出 器 位 移 在 
5~7 mm 之 间 变 化 , 在 不 同 的 排出 器 位 移 下 , 优化 小 
孔 与 双向 进 气 阀 C 评 度 ， 获得 最 低 制冷 温度 


图 9~1 示 , 排出 器 位 移 从 7 mm 至 5 mm 变 

化 ， 最低 制 次 温度 升 高 ， 热 压缩 机 压力 波动 幅 值 减 

小 ,竹山 到 制冷 机 部 分 的 声 功 减 小 ,所 需 预 冷 量 从 

A 减 小 到 17.88 W。 排出 屁 位 移 减 小 , 扫 气 容 
必然 会 产生 上 述 结果 . 


6.0 6.5 7.0 
排出 器 位 移 /mm 
图 9 最 低 制 冷 温 度 随 排出 器 位 移 的 变化 


Fig. 9 Dependency of the lowest no-load temperature on the 
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| -一 压力 波动 幅 值 0.36 

i | & 本 文 基于 Sage 对 VM 型 脉冲 管制 冷 机 进行 了 

| "7 三。 整 机 计算 , 考察 其 在 不 同 运行 频率 、 平均 压力 、 排 出 

"了 ?是 。 器 位 移 下 无 负荷 制冷 温度 、 压 力 波动 幅 值 、 声 切 等 的 

| 030 党 变化 规律 . 在 最 优 的 小 孔 与 双向 进 气 立 开 度 下 , 无 负 

028 站 。 和 蓓 制冷 温度 随 频率 的 增 大 先 降低 后 升 高 ， 存 在 一 个 

1 1 最 佳 频率 ， 随 平均 压力 的 增 大 而 降低 ， 随 排出 器 位 

移 的 增 大 而 降低 .在 改造 后 的 实验 装置 上 进行 了 初 

排 轴 器 位 移 /mm 步 的 实验 验证 , 在 平均 压力 2.2 MPa, 频率 3 Hz, 排 

a 0 出 器 位 移 6 mm, 热 压缩 机 热 端 温度 300 K, 冷 端 温 

图 10 A 冷 端 声 功 。。 言 90 K 的 运行 工 况 下 ， 获得 了 14.4 K 的 无 负荷 制 


Fig. 10 Dependency of pressure wave amplitude and acoustic 
power at warm and cold end for low stage regenerator on the 
displacer displacement amplitude 


3 实验 研究 


本 文 在 改造 热 压 缩 机 回 热 器 后 进行 了 初步 的 驶 入 


机 实验 研究 . 温度 计 分 布 如 图 1 所 示 . T1~ 泌 
冷 用 制冷 机 冷 关 和 热 桥 各 点 温度 (Tl 和 Xe 和 
T8 分 别 对 称 布置 ), T9 为 低温 级 冷 头 六 0 为 低 
温 级 脉 管 冷 端 温 度 ，T11、T12 分 别 葡 肚 管 1/3, 2/3 
处 温度 , 热 压缩 机 运行 参数 如 示 ， 
表 3 热 压 缩 机 运行 参数 
Table 3 The operating parameters of thermal 
compressor 


f/Hz Pn /MPa AP/MPa 
3 2.2 0.21 


Pr 
1.22 


温度 T2, T7( 即 热 压 缩 机 冷 端 温度 ) 达到 70 KC 
时 开启 热 压缩 机 ， 低 漫 级 制冷 机 开始 工作 .需要 指 
出 的 是 , 实验 中 发 现 当 排出 恬 位 移 超过 6 mm 时 , 热 
压缩 机 有 异 响 , 因此 实验 中 最 大 位 移 均 控制 在 6 mm 
以 内 . 通过 不 断 调节 小 孔 、 双 向 进 气 立 门 的 开 度 , 最 
终 获得 了 14.4 KK 的 无 负 蓓 制冷 温度 , 预 冷 用 制冷 机 
提供 的 预 冷 量 大 约 为 49.1W@63.7K。 由 于 实验 过 程 
中 小 孔 、 双 向 进 气 阀 门 的 开 度 需要 不 断 的 优化 调节 ， 
目前 只 进行 了 初步 的 实验 研究 ， 其 它 工 况 的 实验 还 
在 继续 进行 当中 ， 


Ce 


实验 过 程 中 ， 先 开启 预 冷 用 压缩 机 ， 待 其 项 ;3 


冷 温 度 ， 与 相同 工 况 下 的 理论 计算 结果 相当 ， 进 一 
步 的 实验 研究 在 进行 当中 。， 
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